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ABSTRACT 

No apparent correlation was found between giant pulses (GPs) and X-ray photons from the Crab 
pulsar during 5.4 hours of simultaneous observations with the Green Bank Telescope at 1.5 GHz 
and Chandra X-Ray Observatory primarily in the energy range 1.5— 4.5keV. During the Crab pulsar 
periods with GPs the X-ray flux in radio emission phase windows does not change more than by 
±10% for main pulse (MP) GPs and ±30% for interpulse (IP) GPs. During giant pulses themselves, 
the X-ray flux does not change more than by two times for MP GPs and 5 times for IP GPs. All 
limits quoted are compatible with 2-cr fluctuations of the X-ray flux around the sets of false GPs with 
random arrival times. The results speak in favor of changes in plasma coherence as the origin of GPs. 
However, the results do not rule out variations in the rate of particle creation if the particles that 
emit coherent radio emission are mostly at the lowest Landau level. 
Subject headings: Crab pulsar, Giant Pulses, Chandra 



1. INTRODUCTION 

Giant pulses are a peculiar form of pulsar radio emis- 
sion, clearly distinct from the regular pulses. Only a 
handful of pulsars have been reported to have GPs and 
neither the nature of GP emission nor the strict definition 
of a giant pulse are yet fully clear. The "classical" giant 
pulses, from the pulsars with high values of magnetic 
field at the light cylinder, are short (ns - fis), bright (up 
to MJy) and obey a power- law energy distribution (see 
lKuzminll2007l and references therein) . These properties 
obviously point to the coherent nature of the emission 
mechanism. 

It is interesting that for the pulsars with "classical" 
GPs and pulsed X-rays or 7-rays, giant pulses come 
from the same pulse phase range as high -energy emis- 
sion (fKnighl HoM [Cuiu^noeiil] [20031). That indi- 
cates that GPs and high-energy photons may originate 
in the same region in the magnetosphere. Without any 
doubt, the fact that GPs are (or are not) accompanied 
by increasing high-energy flux is important for testing 
the theories of GP generation. On a broad scale, if the 
GP emission mechanism involves some change in parti- 
cle creation rate or beam direction alteration, one should 
observe a correlation between GPs and high-energy emis- 
sion. If GPs are due to changes in coherence in magne- 
tosphcric plasma, there should be no correlation, since 
high-energy emission is of incoherent origin. Also, while 
propagating in the magnetosphere, radio photons can be 
shifted to a higher energy band by interacting with the 
relativistic plasma via inverse Compton scattering or res- 
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onant absorption. This would result in potentially de- 
tectable radio/high-energy correlation. 

The Crab pulsar is the most studied pulsar with GPs 
and the easiest object for such correlation studies, since 
the pulsar is relatively bright at all wavelengths. GPs 
from the Crab appear in phase regions of the main pulse 
(MP) and interpulse (IP) at all radio frequencies where 
the Crab has been observed (from 23 MHz up to 15 GHz). 
Also, the detection of GPs in high f requency component s 
at 8.35 GHz has been reported by IJessner et al.1 (12005ft. 
albei t not confirmed by later observations (|jessner et al.1 
[2010h . 

Around 4 GHz the interpulse disappears and 
then reappears shifted ahead by about 10 degrees 
(jMoffett fc Hankinsl Il99fift . Curiously, above 4 GHz IP 
CPs have totally different spec tral properties than MP 
GPs at the same frequency (jHankins fe Eilekl 12007ft . 
Thus, IP GPs above 4 GHz may be produced by a 
separate emission mechanism. 

Previously, there have been several correlation stud- 
ies between Crab GPs and its high-energy emission 
(see Table [T] for summary). The only positive re- 
sults so far were obtained in the optical light by 
IShearer et al.l (12003ft . They found that optical pulses 
coincident with GPs were on average 3% brighter than 
the others. A number of 7-ray studies covering the 
energy range from 50 McV to 300 GeV did not find 
any statistically significant correlation, showing that 
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Mickaliger et al., in prep). 

A few studies have also been conducted of Crab X-ray 
variability based on high-energy data alone. The analysis 
of 1 hour of RXTE data (1-10 kcV) showed that the 
flux from the Cr ab pulsar is stable at the level of 7% 
(jPatt et al.lll999ft . Also, i n the 13.3 - 58.4 k eV, based 
on 18 hours of RXTE data. iVivekanandl (|2001ft found the 
rms of X-ray flux variation in time at any phase in the on- 
pulse window to be within of «3.5% of the corresponding 
mean flux of the integrated profile of Crab pulsar. 



2 A. V. Bilous et al. 



Table 1 

Previous studies of correlation between GPs from the Crab pulsar and its high-energy emission 
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Figure 1. An example of a strong GP, dedispersed with the 
DM = 56.8065 pc cm" 3 from the Jodrell Bank Crab Pulsar 
Monthly Ephemeris from Nov 15, 2010 (left) and with the cor- 
rect DM, corrected by 0.0102 pc cm -3 (right). The gap between 

1.3 and 1.5 GHz is due to a backend malfunction. 

These studies ruled out the existence of short strong 
"GP-like" bursts of X-ray emission, but they were not 
capable of detecting small sporadic X-ray flux variations, 
as could be expected around GPs. In this paper we used 

5.4 hours of simultaneous observations of the Crab pulsar 
with the 100-m Robert C. Byrd Green Bank Telescope 
(GBT) and the Chandra X-ray Observatory to perform 
an extensive correlation of GPs with X-ray photons. 

2. RADIO OBSERVATIONS 

Radio observations were carried out on 2010 November 
16 with the GBT using the Green Bank Ultimate Pulsar 
Processing Instrument (GUPPI) at a central frequency of 

1.5 GHz in coherent dedispersion search mode. The total 
nominal band of 800 MHz was split into 128 frequency 
channels and the full Stokes parameters were recorded 
with a sampling interval of i rcs = 2.56 fj,s. Unfortunately, 
the data from 1.3 to 1.5 GHz were lost because of a back- 
end malfunction, thus reducing the effective bandwidth 
to 600 MHz. Radio observing time comprised 5.4 hours 
or 5.8 x 10 5 pulsar periods. 



At 1.5 GHz the system temperatu re is dominat ed by 
the Crab nebula, which is smaller (Green 2009) than 
the 8.2' GBT beam at this frequency. We calculated 
the flux density for the Crab Ne bula with the rela - 
tion S(f) = 955 x (//GHz)" - 27 Jy (jCordes et al.ll200l . 
With the receiver temperature T rec = 20 K, background 
sky temperature Tbkg = 1.5 K, and gain G = 2 K/ Jy, our 
estimated system equivalent flux density SEFD = 7.9 Jy 
for a 600 MHz bandwidth and a sampling time of 2.56 ^s. 

During the observations, we performed coherent dedis- 
persion in each of the frequency channels with the most 
up-to-date value of dispersion measure (DM) from the 
Jodrell Bank Crab pulsar monthly ephcmcrifQ, DMjb = 
56.8065 pc cm -3 . However, the dispersion measure of the 
Crab pulsar is known to vary on monthly timescales by 
~ 0.01 pc cm -3 . Such error in DM will not result in 
detectable pulse smearing within one frequency chan- 
nel, but across 800 MHz bandwidth the net smearing 
can reach 20 /is, or 8 time samples. To check the DM, 
we looked at the spectra of the brightest events dedis- 
persed with DMjb. All events exhibited the character- 
istic quadratic trend across the band, as expected for 
an error in DM (see Fig. [TJ left). We dedispersed these 
events with a set of trial DM values and for each pulse 
found the DM which corresponded to maximum peak 
S/N. The mean of these DM values was 56.7963 peem -3 
with standard deviation of 0.001 pc cm -3 . An example 
of a giant pulse dedispersed with the corrected dispersion 
measure is shown in Fig. [TJ (right). 

After correcting the dispersion measure, we dedis- 
persed the raw data using the PRESTO packagfl For 
the correlation study, we used only the total intensity 
time series. More detailed investigation of the recorded 
GPs, including spectra, polarization, width and energy 
distribution will be presented elsewhere. 

To search for GPs we convolved the dedispersed time 
series with a series of boxcar functions of varying boxcar 
width. We used singlepulse_search.py from PRESTO, 
which accounts for all the possible "phase" offsets of each 
boxcar and records all events with S/N>5 on the aver- 

7 http : /Vwww.i b .man. ac .uk/$\sim$pulsar/crab .html see also 
iLvne et alJ (liWJ) 

8 |http: //www. cv.nrao . e du/$\sim$sransom/presto/ 



Correlation of Chandra photons with the radio GPs from the Crab pulsar 



3 




10- 



MP phase range 
IP phase range 

ilse 




0.4 0,6 
Phase 



10' 
Energy, Jy.fjs 



10" 



Figure 2. Left: GP candidates, selected by matched-filtering with a set of boxcars of different width. Most of the candidates come from 
phase regions coincident with the MP (red triangles) and IP (blue crosses). Events outside the MP/IP (grey dots) tend to concentrate at a 
set of levels because of the width-dependent selection thresholds (see text for explanation) . The average profile from the same observation 
is plotted on top for reference. Right: Energy distributions of selected GPs and spurious detections. The energy distribution of spurious 
detections was scaled by 0.023 to match the size of MP/IP phase windows. The shape of the distributions below 1011 Jy-fis is affected by 
width-dependent selection thresholds. 



aged signal. Such definition of the selection threshold fa- 
cilitates sifting the duplicate candidates (the ones above 
the threshold for different boxcar widths, but with the 
same time of arrival). While comparing different box- 
car widths, the peak S/N of the averaged signal has a 
clear maximum when the boxcar width matches the pulse 
width. However, selecting a candidate by its S/N on the 
averaged signal effectively results in a width-dependent 
energy threshold. After averaging by n samples, the stan- 
dard deviation of the signal drops by ^Jn. Therefore S /N 
of 5 on the averaged signal corresponds to a mean S/N of 
5/y/n on the signal with original time resolution. The en- 
ergy of each candidate was calculated by multiplying the 
mean S/N by the corresponding boxcar width. Thus, the 
minimum detectable energy of the n-sample wide candi- 
date was 5\fn ■ SEFD • i rcs , or 101. ly^ Jy'M s wrtn SEFD 
of 7.9 Jy and time resolution of 2.56 /is. In our analysis 
we used boxcars of 1, 2, 3, 4, 6, 9, 14, 20, 30, 45, 70 and 
100 samples and the width-dependent energy thresholds 
can be seen on Fig. [2] (left) as the concentration of events 
at the energy levels proportional to the square root of the 
boxcar widths used. 

Most of GP candidates were detected within two nar- 
row phase regions, each spanning 0.022 pulsar rotations 
and coincident with the main pulse and the intcrpulse on 
the average profile, obtained by folding the dedispersed 
time series with the pulsar period (see Fig.[5J left). Out- 
side the MP/IP region the phase distribution of events 
is flat, as it should be for noise. No excess of events was 
found in the Low Frequency Component (peaks at the 
phase 0.05). 

As has been reported bv lPopov fc Stappersl (j2007t ). gi- 
ant pulses from the MP and IP have different energy 
distributions: IP GPs are well fitted by a single power- 
law, whereas MP GPs obey a broken power-law, with the 
energy of the break depending on observing frequency as 
-^break = 7(^/GHz) 3 ' 4 kJy-^s. For our observing fre- 
quencies, inbreak ranges from 0.8 to 5.0kJy-/xs, making 
the distribution of MP GPs detected in the 800 MHz 



band a convolution of several different broken power- 
laws. 

As can be seen from Fig. [5] (right), above 300 Jy-^s the 
distributions for events in MP and IP windows behave as 
predicted. We believe that the roll-off below 300 Jy-^s is 
most probably due to the underestimation of width (and 
thus energy) for faint pulses. Also, all distributions are 
affected by width-dependent selection thresholds. This 
impact is best visible as a dip between 101.1Jy-/US and 
lOl.lv^Jy-Ats. In any case, the lowest energy of events is 
20 (for MP window) or 100 (for IP) times bigger than the 
corresponding energy of the mean profile. By the formal 
definition of GPs as those with e nergy exceedin g 10 times 
the energy of the mean profile (Cairns 2004), even the 
faintest pulses in our sample can be considered giant. 
Thus, for the subsequent analysis we kept all events in 
the on-pulse phase range, 25986 in MP and 4636 in IP. 
In the off-pulse phase range, which spanned 0.956 pulsar 
rotations, 23366 events were detected. Thus, we estimate 
the fraction of false GPs to be about 2.1% in the MP 
window and 11.6% in the IP window. 

W e put the list of G P arrival times in TEMP02 for- 
mat dHobbs et al.ll2006l ) and converted to the barycentric 
reference frame for the subsequent correlation analysis 
with Chandra data. We used the TDB time system and 
DE405 solar system ephemeris. Times of arrival (TOAs) 
were corrected for delay due to propagation in the ionized 
interstellar medium (ISM). Timing errors due to errors 
in DM are up to 1.8 /xs. 

The observed rate and peak flux densities of GPs are 
inevitably affected by scintillation in the ISM and the 
Crab Nebula. To find the scintillation timescales we con- 
structed the autocorrelation function for the GP rate. 
We did not find any GP rate change on any timescales 
down to 1 s, the smallest bin size possible with our lim- 
ited sample of GPs. Also, we found no scintillation pat- 
terns in the two-dimensional (observing frequency and 
time) autocor relation function for the dedispersed data 
(IRickettlll97"l . averaged over 5 seconds to increase the 
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Figure 3. Left: central part of the HRC-S image. Solid line: source region, zero order (circle) and dispersed spectra (boxes, only central 
part shown). Dashed line: background region. Due to an imperfect PSF model, we still detect weak pulsed emission in the background, so 
we used only photons within pulsar phases 0.7—1.0. Right: phase distribution of the extracted X-ray photons (black). The mean profile of 
the radio emission is overplotted for reference (grey). 



signal to noise ratio. 

By fitting an exponential scattering tail to the average 
profile of the strongest pulses (S/N> 40), we estimated 
the scattering time to be about 10.3 fj,s, which is about 
10 times bigger t han the expected br oadening by inter- 
stellar scattering ([Kuzmin et al. 2002]) and could be ex- 
plained by scintil lations of the nebula or by the intrinsic 
duration of GPs (|Crosslev et al.ll2010Q . In any case, both 
expected and measured scattering times give a decorre- 
lation bandwidth at least an order of magnitude smaller 
than the width of one frequency channel. The refractive 
interstellar scintillation (HISS) timcscale, triss, can be 
scaled using the v ~ 2 2 dependence deriv ed from a five 
frequency dataset (|Rickett fc Lvnel fl990i) . For 1.5 GHz 
it corresponds to Triss=1-2 days, much bigger than our 
observing timcspan. The diffractive interstellar timescale 
is tdiss = triss • A /diss// = Is. However, we did not 
have enough GPs or sensitivity in folded data series to 
detect scintillations on such short timescales. 

3. X-RAY OBSERVATIONS 

The X-ray observations were carried out with Chandra, 
using the Low-Energy Transmission Grating (LETG) 
and High-Resolution Camera spectroscopy detector 
(HRC-S) in timing mode (Obsld 11245). To ensure good 
timing accuracy we used a specific HRC set-up that se- 
lected only signals from the central HRC-S segment, dis- 
abled all on-board vetoing, and set the trigger-level = 
50. The errors in photon arrival times in timing mode 
are about 16 Simultaneous observing time with the 
GBT comprised 5.4 hours. 

We created the level=2 photon lists using the data from 
the Chandra X-ray Center pipelinj TO l. and barycentered 
the arrival times with axbary from the CIAO software 
package. Axbary converts all arrival times to the TDB 
time system, so we could assign each photon its pulse 
phase using TEMP02 with the same pulsar and solar 
system ephemeris as for the radio data. 



9 |http: //asc .harvard.edu/cal/Hrc/ 1 
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While checking the X-ray data, we noticed that the 
count flux from any region on the detector gradually 
grows by about 20% over the course of observations. 
Specifically, if f(t) is the count flux in some region at a 

certain moment, then f(t) can be expressed as k(t) ■ f(t), 
where f(t) is the average count flux in that region and 
k(t) has roughly the same shape for all regions. Inter- 
estingly, the phase-resolved flux exhibits the same be- 
havior - the increase of flux in each phase bin is pro- 
portional to the mean value of flux in that phase bin. 
Such a variability pattern cannot be explained by in- 
creasing background, since its change would be indepen- 
dent from the pulsar phase. Discussions with the HRC 
instrument engineers and scientists did not result in an 
explanation for the rate change (CXC Help Desk Ticket 
#13509). The detector particle background appeared to 
be essentially flat throughout the observation. The only 
change in instrument conditions seen was a slight tem- 
perature change. However, we have no mechanism to 
explain how it could change the X-ray rate but not the 
charged-particle rate (Michael Juda, private communica- 
tion). 

Since X-ray flux changes gradually and GP rate is con- 
stant over the course of observations, we can still con- 
strain the change in X-ray flux around GPs, by compar- 
ing the X-ray flux around GPs to the X-ray flux around 
random time points with the same pulsar spin phase 
range as for the GPs. However, measuring the absolute 
value of the flux is currently not possible. 

For the subsequent analysis we extracted photons from 
1.9" radius circle around the pulsar image in the 0th 
order LETG spectrum and from the two boxes around 
higher spectral orders (see Fig. [31 left, solid line). We 
did not make any cuts by energy, but most photons in 
our sample had energies of 1.5— 4.5keV. Figure 02 right 
shows the phase distribution for our X-ray sample, with 
the radio profile overplotted for reference. Giant pulses 
(which coincide with the peak of radio emission) are de- 
layed with respect to the maximum of the X-ray profile 
by about 0. 01 of pulsar phase . This confirms the result 
reported bv lRots et all (|2004f ). 
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Figure 4. Upper subplots: X-ray counts in the radio emission region around MP GPs (left, within dashed lines) and IP GPs (right, within 
dashed lines) together with X-ray counts in the windows of similar size, shifted by the size of the window within one pulsar period from a 
GP. Black line: X-ray counts around the same number of MP/IP GPs with randomized integer part of phase. Shaded regions - 68% and 
95% percentile on the pool of simulated datasets. Bottom subplots: relative change of X-ray flux (/gp — /noGP )//noGP ■ Only regions with 
good statistics are shown. 



The average count rate on the whole detector was 
about 100 cts/s, which is well below telemetry satura- 
tion limit of 184 cts/s. The average dead time fraction 
was 0.9897, so the dead time after each recorded photon 
comprised 100 fis. This time is bigger than the duration 
of a typical GP and only a few times smaller than the size 
of MP/IP emission window. However, since the average 
number of X-ray photons per one pulsar period in our ob- 
servations is quite low - only 3 photons from the whole 
detector and 0.2 photons from our region of interest, the 
dead time will not considerably affect the correlation re- 
sults on any time scale, unless the correlated photons 
come in bunches of less than 100 /xs. 

4. CORRELATION ANALYSIS 

Ideally, for measuring the change in X-ray flux in the 
vicinity of GPs one should identify all GPs in the data 
sample and compare the X-ray flux in a selected window 
close to GPs (/gp) to the X-ray flux in an identical phase 
window close to the periods without GPs (/ no Gp)- 

In practice, our GP samples always contain some frac- 
tion of spurious detections (s). We do not know the true 
fraction of GPs in the sample (£) and there is background 
X-ray flux (/bkg)- 

Thus, if we measure the X-ray flux in some phase win- 
dow close to GPs, the measured flux will be related to 
the true flux /gp as follows (assuming that s and t are 
small) : 



/micas 
/GP 



/bkg + (1 - s)/gp + s/ 



noGP- 



:i) 



At the same time, by counting photons in the same 
window close to the randomly selected periods we can 
estimate the average (over all periods) X-ray flux in our 
correlation window: 



fwcT = /bkg + i/GP + (1 - 0/noGP- 



(2) 



Then, /gp and / no GP can be found by solving both 
equations simultaneously. 



We performed this procedure separately for MP and 
IP GPs and for the three different timescales - GP du- 
ration, radio emission phase window and pulsar period. 
For every timescale we used not only the correlation win- 
dow centered on GPs, but also a set of windows of the 
same size, shifted with respect to GPs. 

Background flux /bkg was measured in the narrow re- 
gion around the initial extraction region (Fig. [3l left, 
dashed line). Due to imperfect PSF modeling, we still 
detected weak pulsed emission in the background, so we 
used only photons with pulse phases 0.7—1.0. The frac- 
tion of false GPs s was calculated as in Sec. [2] For the 
fraction of periods with GPs, we adopted the fraction of 
periods with detected GPs, 0.04 for MP GPs and 0.008 
for IP GPs, although the true fraction remains unknown. 
However, the limits on flux change are not sensitive to 
the change of t, as long as it remains small. For the two 
larger correlation timescales, the radio emission phase 
window and the pulsar period, varying t up to 0.2 does 
not increase the estimated X-ray flux change in vicinity 
of GP by more than a factor of two. 

For all correlation windows, for each GP we selected 
all X-ray photons within the given window (photons that 
satisfied the criterion for multiple GPs were counted mul- 
tiple times). For comparison, we also measured f^r S m 
the identical windows by creating 1000 lists of fake GPs 
which had the same phase of arrival as real GPs, but 
random arrival period. We used simulated lists instead 
of averaging directly over all periods in order to get the 
errors on the X-ray flux estimates due to the limited num- 
ber of photons/GPs. The histograms of photons around 
real GPs were compared with mean, 68% and 95% per- 
centiles on the pool of simulated datasets and the rel- 
ative flux change, (/gp — /noGp)//noGP was calculated 
with equations [T] and [2j 

For all unshifted correlation windows the value of the 
X-ray flux stayed within the corresponding 95% percentile 
measured on the pool of simulated GP lists with random 
arrival period. In other words, the measured change of 
the X-ray flux in a window around a GP was within 2-a 
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Figure 5. Left: Relative change of X-ray flux (/gp — / n oGp)//noGP m a period around a GP (within dashed lines) and in one-period 
windows shifted by the size of window up to ±15 periods around GPs. Black dots: change of the flux around real GPs. Shaded regions - 
68% and 95% percentile on the pool of GP lists with random arrival times. Right: Relative change of the X-ray flux around a GP on the 
timescale of the GP duration (within dashed lines) and on the same timescale but with phase lags from —0.015 to 0.004. For the average 
emission profile, the peak of the X-ray main pulse is shifted with respect to the peak of average radio pulse by —0.01. 



intrinsic scatter of this value due to limited number of 
GPs/photons. In a few shifted windows the X-ray flux 
exceeded 2-a level, but it always stayed within 3a. 

4.1. Correlation on the timescale of radio emission 
window 

Figure H] (top part of each subplot) shows the X-ray 
counts in the radio emission region around MP GPs (left, 
within dashed lines) and IP GPs (right, within dashed 
lines) together with X-ray counts in the windows of same 
size, shifted by the size of the window within one pulsar 
period from a GP. The bottom part of each subplot shows 
the relative flux change, (/ G p - /noGp)//noGP- The X- 
ray flux in MP/IP emission regions during periods with 
GPs stays within 2a intrinsic scatter of this value due to 
limited number of GPs/photons, which is ±10% for MP 
and ±30% for IP. 

4.2. Correlation on the timescale of pulsar period 

If the process that generates GPs has somewhat less 
prominent long-term influence on X-ray flux, it would 
be easier to detect it with larger correlation windows. 
Figure [5] (left) shows the relative X-ray flux change for 
a timescale of one pulsar period around GPs (within 
dashed lines) and in one-period windows shifted by the 
size of window up to ±15 pulsar periods around GPs. 
The X-ray flux in the same period as GPs stays within 
2-a intrinsic scatter of this value due to limited number 
of GPs/photons, which is ±5% for MP and ±10% for IP. 

4.3. Correlation on the timescale of GP duration 

If X-ray photons are emitted simultaneously with 
GPs, then the change of X-ray flux will reveal itself on 
timescales comparable to the widths of GP. However, 
if the paths of X-rays and radio emission are different 
throughout the magnetosphere, then the X-ray "pulses" 
will be shifted with respect to the radio GPs. Since the 
radio peak is delayed in phase from the X-ray one by 
0.01, we looked for the correlation on the timescale of 



GPs with 20 phase offsets starting from —0.015 with a 
step size of 0.001. 

Most of our GPs have widths less than the X-ray pho- 
ton timing precision (16 fis), so only rough estimates of 
the change of X-ray flux was possible. We looked for 
the correlation in a window around each GP with the 
size of the window equal to the width of that GP plus 
32 yus. For such choice of correlation window s and t de- 
pend in a complex way on the distribution of widths of 
the GPs. However, they bring only minor corrections, so 
we set both of them to 0. The change of the X-ray flux 
(/gp — /noGp)//noGP was estimated from the following 
system of equations: 



/GP 



/bkg 



/GP ' W GP , /noGP • 32/iS 



32/is + wgp 32/is + wgp 



mcas 
aver 



/bkg + /noGP • 



(3) 
(4) 



From Fig. [5] (right) shows that the relative X-ray flux 
on a timescale of the GP duration stays within ±3cr for 
all phase lags chosen. The X-ray flux during GPs (zero 
phase lag) is within ±2er and does not increase more than 
2 times for MP GPs and 5 times for IP GPs. 

5. DISCUSSION 

The lack of significant change of X-ray flux during giant 
pulses should be accounted for by any theory of GP gen- 
eration. Currently, there exist a number of giant pulse 
theories, which place GPs in different parts of the magne- 
tosphere and draw in various emission mechanisms. Un- 
fortunately, none of the theories that we are aware of 
quantitatively describe possible X-ray byproducts of ra- 
dio GPs. Thus, below we will give only some qualitative 
c onsiderations. 

IGil fc Melikidzd (|2004f ) argue that GPs are due to co- 
herent curvature emission of relativistic charged solitons, 
which fly along magnetic field lines in the polar region 
deep inside the light cylinder. The lack of strong syn- 
chrotron X-ray emission accompanying radio GPs sug- 
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gests that the most particles in solitons have zero pitch 
angle with respe ct to magne t ic field . 

The model of IWeatheraill (|1998t) explains GPs using 
spatial collapse of wave packets in strongly turbulent po- 
lar cap plasma. This model describes direct conversion of 
plasma wave turbulence into radio emission, and since it 
does not imply a significant increase of particle density in 
the emitting volume, little correlation with high-energy 
i ncoherent photons is e xpected. 

lHarding et al.l (|2008f ) directly predicts time and phase 
correlation of photons below 200 McV with radio pulses. 
Radio photons undergo resonant absorption by the sec- 
ondary electron-positron plasma on the field lines inte- 
rior to the slot gap. Because of this, charged particles 
are able to maintain significant pitch angles, and emit 
synchroton X-ray photons. However, quantitative limits 
on the energy of the photons and amount of X-ray flux 
change increase are not given. 

Our results rule out strong correlation between radio 
GPs at 1.1—1.9 GHz and X-ray photons of 1.5-4.5keV. 
However, there can still exist weak correlation on a few 
percen t level, similar to t hose observed in the optical 
band (|Shearer et al.l l2003h and/or correlation between 
GPs at some particular radio wavelength and specific 
high-energy band. 
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